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Abstract 



VO \ Abundances and mj_-spectra of strange and other hadronic particles emanating from central 158- 

200 A GeV reactions between nuclei are found to originate from a thermally equilibrated, deconfined 
source in chemical non-equilibrium. Physical freeze-out parameters and physical conditions (pressure, 
specific energy, entropy, and strangeness) are evaluated. Five properties of the source we obtain are as 

l/^ ' expected for direct hadron emission (hadronization) from a deconfined quark-gluon plasma phase. 
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1 Introduction 

^B ■ 

Quark-gluon plasma (QGP) is, by the meaning of these words, a thermally equilibrated state consisting 
JL . of mobile, color charged quarks and gluons. Thermal equilibrium is established faster than the particle 
abundance (chemical) equilibrium and thus in general quark and gluon abundances in QGP can differ from 
*^h their equilibrium Stefan-Boltzmann limit. This in turn impacts the hadronic particle production yields, and 
as we shall see, the chemical non-equilibrium is a key ingredient in the successful data analysis of experiments 
performed at the CERN-SPS in the past decade. We address here results of 200 A GeV Sulphur (S) beam 
interactions with laboratory stationary targets, such as Gold (Au), Tungsten (W) or Lead (Pb) nuclei and, 
Pb-Pb collisions at 158 A GeV. In these interactions the nominal center of momentum (CM) available 
energy is 8.6-9.2 GeV per participating nucleon. 

Considerable refinement of data analysis has occurred since last comprehensive review of the held [|J . 
Our present work includes in particular the following: 

1. We considered aside of strange also light quark (q = u,d) chemical non-equilibrium abundance || 
and introduce along with the statistical strangeness non-equilibrium parameter j s , its light quark 
analogue j q . 

2. Coulomb distortion of the strange quark phase space has been understood ||. 

Incorporating these developments, we accurately describe abundances of strange as well as non-strange 
hadrons, both at central rapidity and in 4-7r-acceptance. We are thus combining in the present analysis 
the strangeness diagnostic tools of dense matter with the notion of entropy enhancement in deconfined 
phase (§. 

As particle emerge from the QGP hadronization, not only their abundance but also their spectral shape 
is of interest. Our analysis considers the impact of explosive radial flow on the spectra of particles at 
high m±. This contributes significant information about the fireball dynamics and the possibly deconfined 
nature of the hadron source. 
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2 The Coulomb effect in the quark fireball 

The diverse statistical chemical parameters that we need to consider will in a self-explanatory way be 
now introduced, considering the concept of strangeness balance: since strangeness can only be made and 
destroyed by hadronic interactions in pairs, the net abundance in the hadronic fireball must vanish. We 
consider a hot gas of free quarks, and evaluate the difference between strange and anti-strange quark 
numbers (net strangeness). The Coulomb potential originating in the initial proton abundance distorts 
slightly the Fermi particle distributions: strange quarks (negative charge) are attracted, while anti-strange 
quarks repelled. Allowing for this slowly changing potential akin to the relativistic Thomas-Fermi phase 
space occupancy model at finite temperature, the energy of a quark depends on both the momentum and 
the Coulomb potential V: 

E p = Jm?+p*-\v. (1) 



It is helpful to see here in first instance the potential V as a square well, Eq. ([!]) within the volume of 
interest. Counting the states in the fireball we obtain: 
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In Eq. (|2|) the subscript Rf on the spatial integral reminds us that only the classically allowed region within 
the fireball is covered in the integration over the level density, g s (= 6) is the quantum degeneracy of strange 
quarks. The magnitude of the charge of strange quarks (Q s /\e\ = 1/3) is shown explicitly, the Coulomb 
potential refers to a negative integer probe charge. 

The fugacity X s of strange quarks enters particles and antiparticles with opposite power, while the 
occupancy parameter 7 S enters both term with same power. For 7 S < 1 one obtains a rather precise result 
for the range of parameters of interest to us (see below) considering the Boltzmann approximation: 
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The meaning of the different factors is now evident. 7 S controls overall abundance of strange quark pairs, 
multiplying the usual Boltzmann thermal factor while A s controls the difference between the number of 
strange and non-strange quarks. Since strangeness is produced as s,s-pair, the value of A s fulfills the 
constraint 
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which is satisfied exactly both in the Boltzmann limit Eq. (^), and for the exact quantum distribution 
Eq. (§), when: 
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Aq is not a fugacity that can be adjusted to satisfy a chemical condition, its value is defined by the applicable 
Coulomb potential V. More generally, in order to account for the Coulomb effect, the quark fugacities within 
a deconfined region should be renamed as follows in order to absorb the Coulomb potential effect: 

A s —* A s = A s Aq , Arf — ► Xd = A^Aq , 



-2/3 x ~\ - J~\ ~\.-\ \-i/6 



A u — > A M = X u Xq , A g — » X q = \j X u Xd — X q Xg . (6) 

Since Qd = Qs = —1/3, the first line is quite evident after the above strangeness discussion, the second 
follows with Q u = +2/3. Note that for a negatively charged strange quark the tilded fugacity Eq. (||) 
contains a factor with positive power 1/3 but the potential that enters the quantity Aq is negative, and 
thus X s < X s . Because Coulomb-effect acts in opposite way on u and d quarks, its net impact on X q is smaller 
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than on X s , and it also acts in the opposite way with X q > X q . To see the relevance of the tilde fugacities for 
light quarks, note that in order to obtain baryon density in QGP one needs to use the tilde-quark fugacity 
to account for the Coulomb potential influence on the phase space. 

It is somewhat unexpected that for the Pb-Pb fireball the Coulomb effect is at all relevant. Recall that 
for a uniform charge distribution within a radius Rf of charge Zf. 
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for r < Rf ; 

(7) 
for r > Rf . 



Choosing Rf = 8 fm, T = 140 MeV, m s = 200 MeV (value of j s is practically irrelevant) we find as solution 
of Eq. (D for (N s - N s ) = for Z { = 150 the value A s = 1.10 (precisely: 1.0983; X s = 1.10 corresponds 
to R{ = 7.87 fm). We will see that both this values within the experimental precision arise from study of 
particle abundances. For the S-W/Pb reactions this Coulomb effect is practically negligible. 

In the past we (and others) have disregarded in the description of the hadronic final state abundances 
the electrical charges and interactions of the produced hadrons. This is a correct way to analyze the 
chemical properties since, as already mentioned, the quantity Aq is not a new fugacity: conservation of 
flavor already assures charge conservation in chemical hadronic reactions, and use of Aj, i = u,d,s exhausts 
all available chemical balance conditions for the abundances of hadronic particles. As shown here, the 
Coulomb deformation of the phase space in the QGP fireball makes it necessary to rethink the implications 
that the final state particle measured fugacities have on the fireball properties. 

3 Freeze-out of hadrons 

The production of hadrons from a QGP fireball occurs mainly by way of quark coalescence and gluon 
fragmentation, and there can be some quark fragmentation as well. We will explicitly consider the recom- 
bination of quarks, but implicitly the gluon fragmentation is accounted for by our allowance for chemical 
nonequilibrium. The relative number of primary particles freezing out from a source is obtained noting 
that the fugacity and phase space occupancy of a composite hadronic particle can be expressed by its 
constituents and that the probability to find all j-components contained within the i-th emitted particle is: 

N i <xl[j j X j e- E ^ T , Xi = U X j, 7i = Hl J - (8) 

jei jei jei 

Experimental data with full phase space coverage, or central rapidity region \y — t/cmI < 0.5, for m± > 
1.5 GeV are considered; recall that the energy of a hadron T is expressed by the spectral parameters m± 
and y as follows, 



E% = zZ E i ' Ei = \l m2 i +P 2 = \J m2 i + P± cosh(2/ - j/cm) , 

where ycM is the center of momentum rapidity of the fireball formed by the colliding nuclei. 

The yield of particles is controlled by the freeze-out temperature Tf. This freeze-out temperature is 
different from the m^-spectral temperature T±, which also comprises the effect of collective matter flows 
originating in the explosive disintegration driven by the internal pressure of compressed hadronic matter. In 
order to model the flow and freeze-out of the fireball surface, one in general needs several new implicit and/or 
explicit parameters. We therefore will make an effort to choose experimental variables (compatible particle 
ratios) which are largely flow independent. This approach also diminishes the influence of heavy resonance 
population — we include in Eq.(||) hadronic states up to M = 2 GeV, and also include quantum statistical 
corrections, allowing for first Bose and Fermi distribution corrections in the phase space content. It is hard 
to check if indeed we succeeded in eliminating the uncertainty about high mass hadron populations. As 
we shall see comparing descriptions which exclude flow with those that include it, our approach is indeed 
largely flow-insensitive. 
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We consider here a simple radial flow model, with freeze-out in CM frame at constant laboratory time, 
implying that causally disconnected domains of the dense matter fireball are synchronized at the instant 
of the collision. Within this approach ||, the spectra and thus also multiplicities of particles emitted are 
obtained replacing the Boltzmann exponential factor in Eq.(g), 

1 r ~tv e • 

6 ~ Ej/T ~" 2^ / d ^ vlv{l + €c ' Pj/Ej)^^^^ 3 '^ » (9) 

where as usual 7„ = 1/^/1 — v£ . Eq. (||) can be intuitively obtained by a Lorentz transformation between an 
observer on the surface of the fireball, and one at rest in the general CM (laboratory) frame. One common 
feature of all flow scenarios is that, at sufficiently high m±, the spectral temperature (inverse slope) T± can 
be derived from the freeze-out temperature Tf with the help of the Doppler formula: 

T ± = Tav(l + v c ) . (10) 

In actual numerical work, we proceed as follows to account for the Doppler effect: for a given pair of 
values Tf and v c , the resulting m± particle spectrum is obtained and analyzed using the spectral shape and 
procedure employed for the particular collision system by the experimental groups, yielding the theoretical 
inverse slope 'temperature' T| . 

Once the parameters Tf, X q , A s , 7 g , 7 S and v c have been determined studying available particle yields, 
and m± slopes, the entire phase space of particles produced is fully characterized within our elaborate 
statistical model. Our model is in fact just an elaboration of the original Fermi model p], in fact all 
we do is to allow hadronic particles to be produced in the manner dictated by the phase space size of 
valance quarks. With the full understanding of the phase space of all hadrons, we can evaluate the physical 
properties of the system at freeze-out, such as, e.g., energy and entropy per baryon, strangeness content. 

4 Results of data analysis 

As noted our analysis requires that we form particle abundance ratios between what we call compatible 
hadrons. We considered for S-Au/W/Pb reactions 18 data points listed in table [l] (of which three comprise 
results with £Ts). For Pb-Pb we address here 15 presently available particle yield ratios listed in table ^ 
(of which four comprise the STs). We believe to have included in our discussion most if not all particle 
multiplicity results available presently. 

The theoretical particle yield results shown columns in tables |l| and y are obtained looking for a set of 
physical parameters which will minimize the difference between theory and experiment. The resulting total 
error of the ratios R is shown at the bottom of tables [y and J2j 

^I^eTJ • (11) 

It is a non-trivial matter to determine the confidence level that goes with the different data analysis 
approaches since some of the results considered are partially redundant, and a few data points can be 
obtained from others by algebraic relations arising in terms of their theoretical definitions; there are two 
types of such relations: 

n + Ti _n i + n/n a _ a I s 

However, due to smallness of the total error found for some of the approaches it is clear without detailed 
analysis that only these are statistically significant. 

In addition to the abundance data, we also explored the transverse mass mj_-spectra when the collective 
flow velocity was allowed in the description, and the bottom line of tables [l] and ^ in columns D„, F v 
includes in these cases the error found in the inverse slope parameter of the spectra. The procedure we 
used is as follows: since within the error the high mj_ strange (anti)baryon inverse slopes are within error, 




Hadron Freeze-out and QGP Hadronization 



Table 1: Particle ratios studied in our analysis for S-W/Pb/Au reactions: experimental results with 
references and kinematic cuts are given, followed by columns showing results for the different strategies of 
analysis B-F. Asterisk * means a predicted result (corresponding data is not fitted). Subscript s implies 
forced strangeness conservation, subscript v implies inclusion of collective flow. The experimental results 
considered are from: 

S. Abatzis et al, WA85 Collaboration, Heavy Ion Physics 4, 79 (1996). 

S.Abatzis et al, WA85 Collaboration, Phys.Lett.B 347, 158 (1995). 

S. Abatzis et al, WA85 Collaboration, Phys.Lett.B 376, 251 (1996). 

I.G.Bearden et al, NA44 Collaboration, Phys. Rev. C 57, 837 (1998). 

D. Rohrich for the NA35 Collaboration, Heavy Ion Physics 4, 71 (1996). 

S-Ag value adopted here: T.Alber et al, NA35 Collaboration, Eur. Phys. J. C 2, 643 (1998). 

A. Iyono et al, EMU05 Collaboration, Nucl.Phys.A 544, 455c (1992) and Y. Takahashi et al, EMU05 

Collaboration, private communication. 



Ratios 


Rcf. 


Cuts [GcV] 


Exp. Data 


B 


C 


D 


D s 


F 


D„ 


F, 


S/A 


1 


1.2 <p± < 3 


0.097 ± 0.006 


0.16 


0.11 


0.099 


0.11 


0.10 


0.11 


0.11 


S/A 


1 


1.2 <p ± < 3 


0.23 ± 0.02 


0.38 


0.23 


0.22 


0.18 


0.22 


0.23 


0.22 


A/A 


1 


1.2 <p±_ < 3 


0.196 ± 0.011 


0.20 


0.20 


0.203 


0.20 


0.20 


0.20 


0.20 


S/S 


1 


1.2 <p ± < 3 


0.47 ± 0.06 


0.48 


0.44 


0.45 


0.33 


0.44 


0.44 


0.43 


fi/fi 


2 


p± > 1.6 


0.57 ± 0.41 


1.18* 


0.96* 


1.01* 


0.55* 


0.98 


1.09* 


1.05* 


n+n 


2 


p_i_ > 1.6 


0.80 ± 0.40 


0.27* 


0.17* 


0.16* 


0.16* 


0.16 


0.13* 


0.13* 


K+/K- 


1 


p_i_ > 0.9 


1.67 ± 0.15 


2.06 


1.78 


1.82 


1.43 


1.80 


1.77 


1.75 


tf°/A 


3 


P± > 1 


1.43 ± 0.10 


1.56 


1.64 


1.41 


1.25 


1.41 


1.38 


1.39 


K°/A 


3 


p_i_ > 1 


6.45 ± 0.61 


7.79 


8.02 


6.96 


6.18 


6.96 


6.81 


6.86 


K°/A 


1 


m± > 1.9 


0.22 ± 0.02 


0.26 


0.28 


0.24 


0.24 


0.24 


0.24 


0.24 


K°/A 


1 


mi > 1.9 


0.87 ± 0.09 


1.30 


1.38 


1.15 


1.20 


1.16 


1.18 


1.17 


S/A 


1 


mj. > 1.9 


0.17 ± 0.01 


0.27 


0.18 


0.17 


0.18 


0.17 


0.16 


0.17 


S/A 


1 


m± > 1.9 


0.38 ± 0.04 


0.64 


0.38 


0.38 


0.30 


0.37 


0.35 


0.35 


n+n 


1 


mi > 2.3 


1.7 ± 0.9 


0.98* 


0.59* 


0.58* 


0.52* 


0.58 


0.72* 


0.75* 


P/P 


4 


Mid-rapidity 


11 ± 2 


11.2 


10.1 


10.6 


7.96 


10.5 


10.6 


10.4 


A/p 


5 


4 7T 


1.2 ± 0.3 


2.50 


1.47 


1.44 


1.15 


1.43 


1.58 


1.66 


p—p 


6 


4 7T 


4.3 ± 0.3 


4.4 


4.2 


4.1 


3.6 


4.1 


4.2 


4.2 


h + -h~ 
h++h~ 


7 


4 7T 


0.124 ± 0.014 


0.11 


0.10 


0.103 


0.09 


0.10 


0.12 


0.12 


xi 


264 


30 


6.5 


38 


12 


6.2 


11 



overlapping we decided to consider just one 'mean' experimental value T± = 235±10 for S-induced reactions 
and T± = 265 ± 15 for Pb-induced reactions. Thus we add one experimental value and one parameter, 
without changing the number of degrees of freedom. Once we find values of If and v c , we evaluate the 
slopes of the theoretical spectra. The resulting theoretical T t J h values are in remarkable agreement with 
experimental Tj_, well beyond what we expected, as is shown in table y. An exception is the fully strange 
O + $7 spectrum. We note in passing that when v c was introduced we found little additional correlation 
between now 6 theoretical parameters. The collective flow velocity is a truly new degree of freedom and it 
helps to attain a more consistent description of the experimental data available. 

Although it is clear that one should be using a full-fledged model such as D v , we address also cases B and 
C. The reason for this arises from our desire to demonstrate the empirical need for chemical non-equilibrium: 
in the approach B, complete chemical equilibrium j± = 1 is assumed. As we see in tables [l] and || this 
approach has rather large error. Despite this the results in column B in tables [l] and || are often compared 
favorably with experiment, indeed this result can be presented quite convincingly on a logarithmic scale. Yet 
as we see the disagreement between theory and experiment is quite forbidding. With this remark we wish to 
demonstrate the need for comprehensive and precisely measured hadron abundance data sample, including 
abundances of multi-strange antibaryons, which were already 20 years ago identified as the hadronic signals 
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Table 2: Particle ratios studied in our analysis for Pb-Pb reactions: experimental results with references 
and kinematic cuts are given, followed by columns showing results for the different strategies of analysis 
B-F. Asterisk * means a predicted result (corresponding data is not fitted or not available). Subscript s 
implies forced strangeness conservation, subscript v implies inclusion of collective flow. The experimental 
results considered are from: 

1 I.Kralik, for the WA97 Collaboration, Nucl. Phys. A 638,115, (1998). 

2 G.J. Odyniec, for the NA49 Collaboration, J. Phys. G 23, 1827 (1997). 

3 P.G.Jones, for the NA49 Collaboration, Nucl. Phys. A 610, 188c (1996). 

4 F.Piihlhofer, for the NA49 Collaboration, Nucl. Phys. A 638, 431,(1998). 

5 C.Bormann, for the NA49 Collaboration, J. Phys. G 23, 1817 (1997). 

6 G.J. Odyniec, Nucl. Phys. A 638, 135, (1998). 

7 D. Rohrig, for the NA49 Collaboration, "Recent results from NA49 experiment on Pb-Pb collisions at 
158 A GeV", see Fig. 4, in proc. of EPS-HEP Conference, Jerusalem, Aug. 19-26, 1997. 

8 A. K. Holme, for the WA97 Collaboration, J. Phys. G 23, 1851 (1997). 



Ratios 


Ref. 


Cuts [GeV] 


Exp. Data 


B 


C 


D 


D s 


F 


D„ 


F« 


S/A 


1 


p± > 0.7 


0.099 ± 0.008 


0.138 


0.093 


0.095 


0.098 


0.107 


0.102 


0.110 


S/A 


1 


p±_ > 0.7 


0.203 ± 0.024 


0.322 


0.198 


0.206 


0.215 


0.216 


0.210 


0.195 


A/A 


1 


p± > 0.7 


0.124 ± 0.013 


0.100 


0.121 


0.120 


0.119 


0.121 


0.123 


0.128 


S/S 


1 


p ± > 0.7 


0.255 ± 0.025 


0.232 


0.258 


0.260 


0.263 


0.246 


0.252 


0.225 


+ + 


2 


P± > 1. 


0.13 ± 0.03 


0.169 


0.114 


0.118 


0.122 


0.120 


0.123 


0.121 


K°J4> 


3.4 




11.9 ± 1.5 


6.3 


10.4 


9.89 


9.69 


16.1 


12.9 


15.1 


K+/K- 


5 




1.80 ± 0.10 


1.96 


1.75 


1.76 


1.73 


1.62 


1.87 


1.56 


P/P 


6 




18.1 ±4. 


22.0 


17.1 


17.3 


17.9 


16.7 


17.4 


15.3 


K/p 


7 




3. ± 1. 


3.02 


2.91 


2.68 


3.45 


0.65 


2.02 


1.29 


K°/B 


3 




0.183 ± 0.027 


0.305 


0.224 


0.194 


0.167 


0.242 


0.201 


0.281 


h-/B 


3 




1.83 ± 0.2 


1.47 


1.59 


1.80 


1.86 


1.27 


1.83 


1.55 


n/s 


1 


p± > 0.7 


0.192 ± 0.024 


0.119* 


0.080* 


0.078* 


0.080* 


0.192 


0.077* 


0.190 


fi/E 


8 


p± > 0.7 


0.27 ± 0.06 


0.28* 


0.17* 


0.17* 


0.18* 


0.40 


0.18* 


0.40 


n/n 


1 


p± > 0.7 


0.38 ± 0.10 


0.55* 


0.56* 


0.57* 


0.59* 


0.51 


0.60* 


0.47 


(S+S) 


8 


p±_ > 0.7 


0.20 ± 0.03 


0.15* 


0.10* 


0.10* 


0.10* 


0.23 


0.10* 


0.23 


x\ 


88 


24 


1.6 


2.7 


19 


1.5 


18 



Table 3: Particle spectra inverse slopes: theoretical values T t h are obtained imitating the experimental 
procedure from the Tf, u c -parameters. Top portion: S-W experimental T± from experiment WA85 for 
kaons, lambdas and cascades; bottom portion: experimental Pb-Pb Tj_ from experiment NA49 for kaons 
and from experiment WA97 for baryons. The experimental results are from: 

D. Evans, for the WA85-collaboration, APH N.S., Heavy Ion Physics 4, 79 (1996). 

E. Andersen et at, WA97-collaboration, Phys. Lett. B 433, 209, (1998). 

S. Margetis, for the NA49-collaboration, J. Physics G, Nucl. and Part. Phys. 25, 189 (1999). 





T K° 


rpk 


T^A 


r s 


T s 


rpn+n 


Tx [MeV] 
T th [MeV] 


219 ± 5 
215 


233 ± 3 
236 


232 ± 7 
236 


244 ± 12 
246 


238 ± 16 
246 


260 


Tx [MeV] 
T th [MeV] 


223 ± 13 

241 


291 ± 18 
280 


280 ± 20 
280 


289 ± 12 
298 


269 ± 22 
298 


237 ± 24 
335 



of QGP phenomena [M]. The strange antibaryon enhancement reported by the experiment WA97 fully 
confirms the role played by these particles ||. 

In the approach C, we introduce strangeness chemical non-equilibrium S, i.e., we also vary j s , keeping 
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Table 4: Statistical parameters which best describe the experimental S-Au/W/Pb results shown in table |l|. 
Asterisk (*) means a fixed (input) value or result of a constraint. In approaches B to D, particle abundance 
ratios comprising 0, are not considered. In case D s strangeness conservation in the particle yields was 
enforced. In case F the three data-points with Q are considered. Lower index v implies that radial collective 
flow velocity has been considered. 



s-w 


T { [MeV] 


\ 


A s 


7s llq 


ll 


V c 


xi 


B 


144 ± 2 


1.53 ± 0.02 


0.97 ± 0.02 


1* 


1* 


0* 


264 


c 


147 ± 2 


1.49 ± 0.02 


1.01 ± 0.02 


0.62 ± 0.02 


1* 


0* 


30 


D 


143 ± 3 


1.50 ± 0.02 


1.00 ± 0.02 


0.60 ± 0.02 


1.22 ± 0.06 


0* 


6.5 


D s 


153 ± 3 


1.42 ± 0.02 


1.10* ± 0.02 


0.56 ± 0.02 


1.26 ± 0.06 


0* 


38 


F 


144 ± 3 


1.49 ± 0.02 


1.00 ± 0.02 


0.60 ± 0.02 


1.22 ± 0.06 


0* 


12 


D„ 


144 ± 2 


1.51 ± 0.02 


1.00 ± 0.02 


0.69 ± 0.03 


1.41 ± 0.08 


0.49 ± 0.02 


6.2 


F« 


145± 2 


1.50 ± 0.02 


0.99 ± 0.02 


0.69 ± 0.03 


1.43 ± 0.08 


0.50 ± 0.02 


11 



Table 5: Statistical parameters which best describe the experimental Pb-Pb results shown in table |2|. 
Asterisk (*) means a fixed (input) value, or result of a constraint. In approaches B to D, particle abundance 
ratios comprising O, are not considered. In case D s strangeness conservation in the particle yields was 
enforced. In case F the four data-points with Q are considered. Lower index v implies that radial flow 
velocity has been considered. 



Pb-Pb 


T{[MeV] 


\ 


A s 


Is Hq 


lq 


V c 


x{ 


B 


142 ± 3 


1.70 ± 0.03 


1.10 ± 0.02 


1* 


1* 


0* 


88 


C 


144 ± 4 


1.62 ± 0.03 


1.10 ± 0.02 


0.63 ± 0.04 


1* 


0* 


24 


D 


134 ± 3 


1.62 ± 0.03 


1.10 ± 0.02 


0.69 ± 0.08 


1.84 ± 0.30 


0* 


1.6 


D s 


133 ± 3 


1.63 ± 0.03 


1.09* ± 0.02 


0.72 ± 0.12 


2.75 ± 0.35 


0* 


2.7 


F 


334 ± 18 


1.61 ± 0.03 


1.12 ± 0.02 


0.50 ± 0.01 


0.18 ± 0.02 


0* 


19 


D„ 


144 ± 2 


1.60 ± 0.02 


1.10 ± 0.02 


0.86 ± 0.03 


1.72 ± 0.08 


0.58 ± 0.02 


1.5 


F„ 


328 ± 17 


1.59 ± 0.03 


1.13 ± 0.02 


0.51 ± 0.03 


0.34 ± 0.13 


0.38 ± 0.31 


18 



7« = 1. A nearly valid experimental data description is now possible, and indeed, when the error bars were 
smaller this was a satisfactory approach adapted in many studies Q. However, only the possibility of light 
quark non-equilibrium in fit D produces a statistically significant data description. 

It is interesting to note that a significant degradation of \\ occurs, especially in the Pb-Pb data, when 
we require in column F that the particle ratios comprising f2 and fi-particles are also described. We thus 
conclude that a large fraction of these particles must be made in processes that are not considered in the 
present model. 

Another notable study case is shown in column D s , with strangeness conservation enforced. Remarkably, 
the S-W data, table 111, do not like this constraint. A possible explanation is that for S-induced reactions, the 
particle abundances are obtained at relatively high p± . Thus only a small fraction of all strange particles is 
actually observed, and therefore the overall strangeness is hard to balance. Similar conclusion results also 
when radial flow is explicitly allowed for, with a significant unbalanced strangeness fractions remaining, 
as we shall discuss below. On the other hand, this constraint is relatively easily satisfied for the Pb-Pb 
collision results, table g, where a much greater proportion of all strange particles is actually experimentally 
detected. 

The statistical parameters associated with the particle abundances described above are shown in the 
table |I] for S-W/Pb and in table [5] for Pb-Pb reactions. The errors of the statistical parameters shown are 
those provided by the program MINUIT96.03 from CERN program library. When the theory describes the 



8 J. Rafelski and J. Letessier 

Table 6: Tf and physical properties (specific energy, entropy, anti-strangeness, net strangeness, pressure 
and volume) of the full hadron phase space characterized by the statistical parameters given in table |I] for 
the reactions S-Au/W/Pb. Asterisk * means fixed input. 



s-w 


T t [MeV] 


Et/B 


St/B 


st/B 


(h - s{) /B 


Pi [GeV/fm 3 ] 


B 


144 ± 2 


8.9 ± 0.5 


50± 3 


1.66 ± 0.06 


0.44 ± 0.02 


0.056 ± 0.005 


c 


147 ± 2 


9.3 ± 0.5 


49± 3 


1.05 ± 0.05 


0.23 ± 0.02 


0.059 ± 0.005 


D 


143 ± 3 


9.1 ± 0.5 


48± 3 


0.91 ± 0.04 


0.20 ± 0.02 


0.082 ± 0.006 


D s 


153 ± 3 


8.9 ± 0.5 


45 ± 3 


0.76 ± 0.04 


0* 


0.133 ± 0.008 


F 


144 ± 2 


9.1 ± 0.5 


48 ± 3 


0.91 ± 0.05 


0.20 ± 0.02 


0.082 ± 0.006 


D« 


144 ± 2 


8.2 ± 0.5 


44 ± 3 


0.72 ± 0.04 


0.18 ± 0.02 


0.124 ± 0.007 


F„ 


145 ± 2 


8.2 ± 0.5 


44 ± 3 


0.73 ± 0.05 


0.17 ± 0.02 


0.123 ± 0.007 



data well, this is a one standard deviation error in theoretical parameters arising from the experimental 
measurement error. In tables Q and [5] in cases in which 7, / 1 we present the ratio 7 S /jq, which corresponds 
approximately to the parameter 7^ in the data studies in which j„ = 1 has been assumed. It is notable that 
whenever we allow phase space occupancy to vary from the equilibrium, a significant deviation is found. In 
the S-W case, table f|, there is a 25% excess in the light quark occupancy, while strange quarks are 25% 
below equilibrium. In Pb-Pb case, table || the ratio of the nonequilibrium parameters 7s/7<? — 0.7 also 
varies little (excluding the failing cases F with fi,il data), though the individual values j s , 7 g can change 
significantly, even between the high confidence cases. 

We note that in the Pb-Pb reaction, table S, 7 S > 1. This is an important finding, since an explanation 
of this effect involves formation prior to freeze-out in the matter at high density of near chemical equilibrium, 
7 S (t < tf) ~ 1. The ongoing rapid expansion (note that the collective velocity at freeze-out is found to 
be l/\/3) preserves this high strangeness yield, and thus we find the result 7 S > 1. In other words the 
strangeness production freeze-out temperature T s > Tf. Thus the strangeness equilibration time is proven 
implicitly to be of magnitude expected in earlier studies of the QGP processes J10|| . It is hard, if not really 
impossible, to arrive at this result without the QGP hypothesis. Moreover, inspecting figure 38 in |n| we 
see that the yield of strangeness we expect from the kinetic theory in QGP is at the level of 0.75 per baryon, 
the level we indeed will determine below. 

Another notable results is A s ~ A s — 1.0 in the S-Au/W/Pb case, see table ^, and A s ~ 1.1 in the Pb-Pb 
case, see table [5], implying here A s = 1, see section ||. We see clearly for both S- and Pb-induced reactions 
a value of X s , characteristic for a source of freely movable strange quarks with balancing strangeness. 



5 Physical properties of the fireball at chemical freeze-out 

Given the precise statistical information about the properties of the hadron phase space, we can determine 
the physical properties of the hadronic particles, see tables g and 0. We show for the same study cases 
B-F, along with their temperature, the specific energy and entropy content per baryon, and specific anti- 
strangeness content, along with specific strangeness asymmetry, and finally pressure at freeze-out. We note 
that it is improper in general to refer to these properties as those of a 'hadronic gas' formed in nuclear 
collisions, as the particles considered may be emitted in sequence from a deconfined source, and thus there 
may never be a evolution stage corresponding to a hadron gas phase. However, the properties presented 
are those carried away by the emitted particles, and thus characterize the properties of their source. 

The energy per baryon seen in the emitted hadrons is, within error, equal to the available specific energy 
of the collision (8.6 GeV for Pb-Pb and about 8.8-8.9 GeV for S-Au/W/Pb). This implies that the fraction 
of energy left in the central fireball must be the same as the fraction of baryon number. We further note 
that hadrons emitted at freeze-out carry away a specific s content which is determined to be 0.72 ± 0.04 
in S-Au/W/Pb case and 0.78 ± 0.04 for the Pb-Pb collisions (cases D v ). Here we see the most significant 
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Table 7: Tf and physical properties (specific energy, entropy, anti-strangeness, net strangeness, pressure 
and volume) of the full hadron phase space characterized by the statistical parameters given in table |5| for 
the reactions Pb-Pb. Asterisk * means fixed input. 



Pb-Pb 


T f [MeV] 


Et/B 


St/B 


at/B 


(If - s{) /B 


Pi [GeV/fm y ] 


B 


142 ± 3 


7.1 ± 0.5 


41 ± 3 


1.02 ± 0.05 


0.21 ± 0.02 


0.053 ± 0.005 


C 


144 ± 4 


7.7 ± 0.5 


42 ± 3 


0.70 ± 0.05 


0.14 ± 0.02 


0.053 ± 0.005 


D 


134 ± 3 


8.3 ± 0.5 


47 ± 3 


0.61 ± 0.04 


0.08 ± 0.01 


0.185 ± 0.012 


D s 


133 ± 3 


8.7 ± 0.5 


48 ± 3 


0.51 ± 0.04 


0* 


0.687 ± 0.030 


F 


334 ± 18 


9.8 ± 0.5 


24 ± 2 


0.78 ± 0.05 


0.06 ± 0.01 


1.64 ± 0.06 


D„ 


144 ± 2 


7.0 ± 0.5 


38 ± 3 


0.78 ± 0.04 


0.01 ± 0.01 


0.247 ± 0.007 


F v 


328 ± 17 


11.2 ± 1.5 


28 ± 3 


0.90 ± 0.05 


0.09 ± 0.02 


1.40 ± 0.06 



impact of flow, as without it the specific strangeness content seemed to diminish as we moved to the larger 
collision system. We have already alluded repeatedly to the fact that for S-Au/W/Pb case the balance of 
strangeness is not seen in the particles observed experimentally. The asymmetry is 18%, with the excess 
emission of s containing hadrons at high p± > 1 GeV. In the Pb-Pb data this effect disappears, perhaps 
since the p± lower cut-off is smaller. One could also imagine that longitudinal flow which is stronger in 
S-Au/W/Pb is responsible for this effect. 

The small reduction of the specific entropy in Pb-Pb compared to the lighter S-Au/W/Pb is driven 
by the greater baryon stopping in the larger system, also seen in the smaller energy per baryon content. 
Both systems freeze out at E/S = 0.185 GeV (energy per unit of entropy). Aside of Tf, this is a second 
universality feature at hadronization of both systems. The overall high specific entropy content agrees well 
with the entropy content evaluation made earlier yj for the S-W case. This is so because the strange 
particle data are indeed fully consistent within the chemical-nonequilibrium description with the 4ir total 
particle multiplicity results. 



6 Current status and conclusions 

We have presented detailed analysis of hadron abundances observed in central S-Au/W/Pb 200 A GeV 
and Pb-Pb 158 A GeV interactions within thermal equilibrium and chemical non-equilibrium phase space 
model of strange and non-strange hadronic particles. In the analysis of the freeze-out the structure of the 
particle source was irrelevant. However, the results that we found for the statistical parameters point rather 
clearly towards a deconfined QGP source, with quark abundances near but not at chemical equilibrium: 
statistical parameters obtained characterize a strange particle source which, both for S-Au/W/Pb and for 
Pb-Pb case, when allowing for Coulomb deformation of the strange and anti-strange quarks, is exactly 
symmetric between s and s quark carriers, as is natural for a deconfined state. 

There are similarities in freeze-out properties seen comparing results presented in tables y and which 
suggest universality in the extensive physical properties of the two freeze-out systems ||. Despite consider- 
ably varying statistical parameters we obtain similar physical conditions such as Tf and E/S corresponding 
possibly to the common physical properties of QGP at its breakup into hadrons. The precision of the 
data description we have reached, see tables [l}-|| strongly suggests that despite its simplicity the model 
we developed to analyze experimental data provides a reliable image of the hadron production processes. 
When flow is allowed for the freeze-out temperature is identical in both physical systems considered, even 
though, e.g., the baryochemical potential jjlb = 3TlnA g and other physical parameters found are slightly 
different. It is worth to restate the values we obtained, case T) v (6 parameter description, no f2, £1): 

T f = 144 ± 2 MeV, fi B = 178 ± 5MeV , for S-Pb/W , 

T f = 144 ± 2 MeV, /z B = 203 ± 5MeV , for Pb-Pb , 
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Even though there is still considerable uncertainty about other freeze-out flow effects, such as longi- 
tudinal flow (memory of the collision axis), the level of consistency and quality of agreement between a 
wide range of experimental data and our chemical non-equilibrium, thermal equilibrium statistical model 
suggests that, for the observables considered, these effects do not matter. 

The key results we found analyzing experimental data are: 

1. X s = 1 for S and Pb collisions ; 

2. 7s Pb >i, 7 ? >i; 

3. uf> = l/v/3 

4. S/B ~ 40 ; 

5. s/5~0.75; 

are in remarkable agreement with the properties of a deconfined QGP source hadronizing without chemical 
reequilibration. The only natural interpretation of our findings is thus that hadronic particles seen at 158- 
200 GeV A nuclear collisions at CERN-SPS are emerging directly from a deconfined QGP state and do not 
undergo a chemical re-equilibration after they have been produced. 
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